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ABSTRACT 

We explore the parameter space of the two temperature pseudo-Newtonian 
Keplerian accretion flow model for the millimeter and shorter wavelength emis- 
sion from Sagittarius A*. A general relativistic ray-tracing code is used to treat 
the radiative transfer of polarized synchrotron emission from the flow. The syn- 
chrotron self-Comptonization and bremsstrahlung emission components are also 
included. It is shown that the model can readily account for the millimeter to 
sub-millimeter emission characteristics with an accretion rate of ~ 6 x 10 17 g ■ s _1 
and an inclination angle of ~ 40°. However, the corresponding model predicted 
near-infrared and X-ray fluxes are more than one order of magnitude lower than 
the observed 'quiescent' state values. While the extended quiescent-state X-ray 
emission has been attributed to thermal emission from the large-scale accretion 
flow, the NIR emission and flares are likely dominated by emission regions either 
within the last stable orbit of a Schwarzschild black hole or associated with out- 
flows. With the viscous parameter derived from numerical simulations, there is 
still a degeneracy between the electron heating rate and the magnetic parameter. 
A fully general relativistic treatment with the black hole spin incorporated will 
resolve these issues. 

Subject headings: black hole physics — Galaxy: center — plasmas — polarization 
— radiative transfer — sub-millimeter 

1 Key Laboratory for Research in Galaxies and Cosmology, The University of Sciences and Technology of 
China, Chinese Academy of Sciences, Hefei 230026, China; mlhuang@ustc.edu.cn, yfyuan@ustc.edu.cn 

2 Department of Physics and Astronomy, University of Glasgow, Glasgow G12 8QQ, UK; 
sliu@astro.gla.ac.uk 

3 Key Laboratory for Research in Galaxies and Cosmology, Shanghai Astronomical Observatory, Chinese 
Academy of Sciences, Shanghai 200030, China; zshen@shao.ac.cn 

4 Academia Sinica, Institute of Astronomy and Astrophysics, Taipei 106, Taiwan; 
mike@asiaa.sinica.edu.tw 

5 Los Alamos National Laboratory, Los Alamos, NM 87545; hli@lanl.gov, fryer@lanl.gov 

6 Physics Department, The University of Arizona, Tucson, AZ 85721 



-2- 



INTRODUCTION 



Sagittarius (Sgr) A*, the compact radio source associated with the super-massive black 
hole at the Galactic center, is perhaps the best source for the study of phy sical processes 
near the event horizon of a black hole (jSchodel et al.l 120021 ; iGhez et al.l 120051 ). Although its 
luminosity is relatively low, with the high resolution and sensitivity of modern instruments, 
the source has been routinely observed from radio to X-rays. It may als o play a role in 



the production of TeV gamma-ray emission from the Galactic center region (ILiu et al.l 12006 



Aharonian et al.ll2006l ). The low luminosity also renders the source optically thin at millime- 
ter and shorter wavelengths. One therefore can observe emission from the very inner region 
close to the black hole directly at these wavelengths. 

Most of Sgr A* emission is released in the sub-millimeter band, where a strong linear 
polarization is observed. The flux density also varies with the variation amplitu de and rate 
decreasing with the decrease of the observation frequency (IHerrnstein et al.l 12004 ) . The vari- 
ation timescale of sub-millimeter emission can be as short as a few hours, slightly longer than 
the dynamical time near the black hole. The source is highly variable at shorter wavelengths 
with the peak luminosity of some near infrared (NIR) and X-ray flares comparable to the 
sub- millimeter luminosity (IBaganoff et al.ll200ll ; iGenzel et al.ll2003l ) The variation timescale 
of these flares is consistent with events occurring within a few Schwarzschild radii of the black 
hole. And correlated flare a ctivities in the X-ray, NIR, millimeter, and radio bands suggest 
outflows during so me flares (jZhao et al.ll2004j ; lYusef-Zadeh et al.ll2007l ; iMarrone et al.l 120071 : 
Eckart et al.ll2008f). The flares also contain rich temporal, sp ectral, and polarization infor- 
mation (IPorquetl 120081 ; iFalanga et al.ll2008l ; lEckart et al.l 120081 ) . Detailed modelling of these 
emission characteristics is expected to probe the geometry near the black hole, processes 
of the general relativistic magnetohydrodynamics (GRMHD) ( IGammie et al.ll20031), and th e 



kinet i cs of e lectron heating and acceleration in a magnetized relativistic plasma (ILiu et al. 
20041 . l2007bh . 



Motiva ted by MHD simulations of radiatively inefficient accretion flows (IHawley fc Balbus 



1 9 9 ll . |2002| ) , we proposed a Keplerian accretio n flow model for the tim e averaged millimeter 
and shorter wavelength emission from Sgr A* (IMelia et al.ll2000l . 120011 ). The model plays im- 
portant roles in co nstraining the time averaged properties of the p lasma near the black hole 



flYuan et al.ll2003f ) for detailed GRMH D studies flNob 



ization characteristics quantitatively (IBromley et al 



e et al.ll2007| N ). It also treats the polar- 



200 ll ). The model is later generali zed 



to a two temperatu re accretion flow in a pseudo-Newtonian potential ( Liu et al. 2007al lbh . 
Huang et al.l (120081 ) first carried out a self-consistent treatment of the transfer of polarized 
emission through the accretion flow and predicted distinct linear and circular polarization 
characteristics between the sub-millimeter and NIR band due to general relativistic light 
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bending and birefringence effects. In this paper, we present the details of the calculations ( 
|2]), explore the model parameter space (§E]), and discuss future developments (§HJ). 



RADIATIVE TRANSFER OF SYNCHROTRON RADIATION 



Although the importance of a self-consistent treatment of polarized ra diative transfe r 



through relativistic plasmas near black holes was recognized a while ago (Me 



there are still significant uncertainties in the quantitative details (jShcherbakov 



rose 



3 ll997h . 



20081 ). The 



thermal synchro tron emission and absorption coefficients have been d erived independently by 
several authors (jLegg fc Westfold 1968 ; Sazonov 1969 ; Melrose 1971 ). which in general show 
agreeme nt. The deriva t ion of the Fa r aday rotation and conversion coefficients was mostly 
done by (lMelrosdll997l ). IShcherbakovl (120081 ) recently derived expressions appropriate for the 
transition from non-relativistic temperatures to relativistic ones. In this section, we have 
a brief discussion of the general theory of polarized radiation transfer and show that the 
Faraday rotation and conversion coefficients are not independent. The radio of the two is 
equal to the ratio of the circular and linear emission coefficients. Therefore our approach 
of deriving the Faraday conversio n coefficient from the emission coefficients and Faraday 
rotation coefficient is appropriate (IHuang et al.ll2008l ). 



2.1. Properties and Transfer of Partially Polarized Radiation 

The synchrotron radiation from an individual particle is elliptically polarized. The major 
axis of the polarization ellipse e 1 is perpendicular to the plane spanned by the magnetic vector 
B and the wave vector k, and the minor axis e 2 oc kxe 1 is perpendicular to the major axis 
and wave vector. The electric field component of the radiation can be projected along the 
major and minor axes, namely the e- and o-component, respectively. Figure Q] shows the 
corresponding geometry in the three dimensional (3D) Cartesian coordinates (x,y,z): 

B/B = (0,1,0), k/k = (0, cos# B , sinfl B ), 

e 1 = (1,0,0), e 2 = (0, sin^ B , -cos^ B ), (1) 

where 9b is the angle between B and k, and k and B are the wave number and the amplitude 
of the magnetic field, respectively. For a population of relativistic particles with a smooth 
pitch angle distribution, the circular polarization (CP) component of synchrotron emission 
almost cancels out, and the emissivity of the e-component is generally much greater than that 
of the o-component. Therefore, the synchrotron radiation from a population of relativistic 
particles is often treated as partially linearly polarized along e 1 . The e-component base e 1 



-4- 



then indicates the electric vector of the linearly polarized (LP) radiation. Since the CP 
component does not cancel out exactly, sometimes it is more convenient to use the right- 
handed and left-handed CP bases 



te 



(2) 



According to the intensity matrix X 1 - 7 defined by iMelrosd ( 119711 ) with an arbitrary or- 
thogonal bases, the total intensity is given by 



trace (I 1 - 



l n + T 



■22 



1 RR +Z 



LL 



(3) 



Then one has the formal polarization vector p 
j2004f ) 



Pu, Pv) introduced by lLandi Degrinnocenti &: Landolf 



Pq 



Pu 



Pv 



J 11 - J 22 
J 11 + 1 22 ' 
J 12 +X 21 
J 11 + X 22 ' 
l (Z l2 -l 21 ] 

J 11 + J 22 



T RR + T LL ' 



(4) 



The total polarization fraction II, CP fraction lie, LP fraction lit,, and the electric vector 
position angle (EVPA) xo are then given by 



n 

tan2xo 



1 - 



4det(J i - 



trace (Jy) 2 

(tf-n 2 .) 1 / 2 : 



1/2 



(P&+P&) 1/2 , 



(5) 



The normalized polarization vector can be rewritten in the form of the 4D Stokes vector as 
S = (I,Q,U,V) T = I(l,pQ,pu,Pv) T '■ The Stokes vector, the intensity matrix, and (J, lie, 
n_L, Xo) ah give a complete description of the properties of partially polarized emission. 



In order to describe the polarized radiative transfer in highly-magnetized plasma, lLandi DegFInnocenti 



( 120041 ) introduce another three 3D formal vector e, if, p, where e is the normalized emission 
vector, which is related to the 4D emission coefficient e = (ej, Eq, eu, £v) T - e = (^q, €u, e v) = 
{sq, Ejj, ey)/I, if = (t)q, t)u, r)y) is the absorption vector, and p is the Faraday rotation vector. 
The total emission and average absorption coefficients are represented by Ej = ejl and rji , 
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respectively. Then the transfer of polarized emission is described with 

dl . _ _ 

fa = -(Vi + V-P-ei) 1 , 

dp _ , _> _ _ _ _ 

— = —f] + vn ■ pip + p x p + e - ejp , 

ds 



(6) 



which govern the evolution of the polarization properties of radiation through a magnetized 
plasma. These equations can be rewritten in the form of Stokes vector as 



ds 



S 



KS, 



(7) 



or explicitly 
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One can use the rotation matrix 



R(x) 



/ l 











V 



cos2x sin2x 
— sin 2\ cos 2% 
1 



(9) 



to transform the Stokes vector from the coordinates (e^e 2 ) to the reference coordinates 
(a, b), where a corresponds to the North at the observer and b corresponds to the East (see 
Figure [T]) . Then we have 



ds 



-S' 



K' S' 



(10) 



where S' = R(x)S, ? = R(x)e, and K' = R(x)KR(~x). 



2.2. Synchrotron Emission, Absorption, and Faraday Coefficients 

For a population of relativistic electrons in the Maxwell distribution, i.e., relativistic 
thermal distribution, with the number density iVo and temperature T e 3> m e c 2 /k, the dis- 
tribution function with respect to the electron energy E is 

N(E) = No 2(W^ eXpi ~ E/kTeh (U) 
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where m e , c, and k are the electron mass, speed of light, and Boltzmann constant , respec- 



1969; 


Pacholczvk 


1970; 


Melrose 


* * 

1971) 
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(12) 



where e is the elementary charge units and xm is the ratio of the emission frequency v to 
the characteristic frequency v c = 3eBsm9 B (kT e ) 2 /47imlc 5 , and the integrations Ii,Iq and 
ly can be calculated with the aid of the modified Bessel functions K a (z): 
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H{x) 
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Iv{xm) 
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(13) 



The emission coefficients along the two axes e 1 and e 2 are given, respectively, by 



El 
E 2 



2 ' 



(14) 



The absorption coefficient can be calculated with the Kirchhoff's law 
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where S = Ei/rjj is the source function, S = B v = 2hu 3 /c 2 [exp (hv/kT e ) — 1] for a thermal 
distribution, and h is the Planck constant. 



According to iMelrosd (119971 ) , the high-frequency waves may be treated as two transverse 
nature wave modes with dispersion relations 

7 2 2 2 1 r 11 , 22 i IV 11 22\2 , ,4 12 21U/21 

k±c — uj = - {a + a ± [(a — a ) + 4a « J ' }, 

where u; = 27T/v, is the wave number, and the response tensor a 1 ^ satisfies a 12 = —a 21 . The 
polarization vectors of the two natural modes are then given by 

± T ± e l +ie 2 ^ a 11 — a 22 =F [(a 11 — a 22 ) 2 + 4a 12 a 21 ] 1 / 2 

e = 7rW ± = 2<^ > ( 16 ) 

which are two orthogonal elliptically polarized modes with axial ratios T±. Note that T + T_ = 
— 1 implying that the two natural modes have identical ellipticity and opposite handedness. 
The electric vector E of any wave can be decomposed on the bases of the two natural modes 
as E + e + + E^e~ or on the bases of the two LP modes as E^e 1 + E^e 2 . Then 

E 1 = E + ^±== + E_ 



VTTt 2 -^T+Tl 

and 

E 2 = E + , = + £_ 



v^Tr 2 _ v / i+^' 

Using the definitions in equation (j2J), a third decomposition can be made on two CP modes 
as E R e R + E L e L , with 

T i + 1 „ T_ + 1 
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The total emission coefficient £/ is proportional to E 2 , i.e., 

+ oc E\ + E 2 _ 
£j = J £i+£ 2 oc E 2 + E 2 . (17) 
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Then the LP and CP emission coefficients (eq and ey) are given by 
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which implies that 



Ey_ 



2ia 



12 



a 11 - a 22 



Prm 
Prrm 



(19) 



where p RM 
Prrm = 
coefficient. 

r] V = e v I S 



la 



ii 



a 



/2ujc is called the rotation measure, or the Faraday rotation coefficient, and 
! )/4u;c is called the relativistic rotation measure, or the Faraday conversion 



Note that e\ = 
= (W+ + W-)/2- 



4± + £2 



V±i 



e±/v± = 5 A Vq = £q/S = (VQ+ +Vq-)/2, 



In a uniform plasma, p — (2p RRM , 0, 2p RM ) and the three formal vector e, ff, and p are 
parallel t o each other. If there is no external radiation, the term (p x 0) in equation s ([6]) 
vanishes (ILandi Degl'Innocenti &: Landolfil 12004 ; lBekefilll966l ; iKennett fc Melroseill998l ) and 
one has the solution with p also parallel to these vectors. Then equations ([6]) become 



dl 
ds" 

dn 

d^ 



-[(i + nn )7 7/ -e / ]/ = - 
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(20) 



where Il r 



+ 4 



+ £ 2 v) 



S I 

i 1 / 2 /e/ is the polarization fraction of an optically thin source. 
Note that in this case, if we denote the unit vector along e as u, then p = Uu, e = Uoeiu, 
fj = el / S = Uorjiu. With the increase of the optical depth, we see that IT decreases, and the 
existence of II enhances the self- absorption by a factor of 1 + nn . 

Therefore, in a uniform plasma, besides the source function that is determined by the 
particle distribution, there are four independent coefficients. From the three emission co- 
efficients given by equations (fT2"l) and the Faraday rotation coefficient py = ^Prm, one can 
derive the absorpti o n and Faraday conversion coefficients with equations f|T5l) and ffT9l) . re- 
spectively. iMelrosd ( 119971 ) derived the three emission coefficients and two Faraday coefficients 
separately. It appears that his results ar e not exactly i n agre ement with equation f|T9l . Nei- 
ther is the p R M and p RR M obtained by IShcherbakovl (120081 ) . All these calculations invoke 
some approximations to simplify the Trubnikov's linear response tensor. While equation 
(|T9|) is derived from very general theoretical considerations without any approximations. We 
therefore stand b y our approach of d eriving p RR M from p RM and the emission coefficients in 
a previous paper (IHuang et al.ll2008l ). 



Melrose! ( 119971 ) gives p R M for thermal plasmas in the cold (7^ = kT e /m e c 2 + 1 ~ 1) and 



extremely relativistic (ER) (7c ^> 1) limits: 

e 3 N B cos 9 B 



Prm 



Prm 



2Ti r y c m 2 c 2 v 2 ' 
In 7 C e 3 N B cos 63 
27 c 2n^ c m 2 c 2 v 2 



(21) 



- 9- 



We use the following extrapolation to obtain /?rm for arbitrary electron temperatures 

Prm = 7" 1 (Prm-Prm)+Prm, (22) 



which is simple and accurate enough according to IShcherbakovl (120081 ). For cold plasmas, 
cyclotron emission dominates, £j > By 3> £q, and py ^> pq. For hot plasmas, synchrotron 
emission dominates, ej > eg ey, and pg 3> py. 



dinates (e 1 



For a given ray, the magnetic field structure of the accretion flow determines the coor- 
e 2 ). For given reference directions of the North and East, we also obtain the 
coordinates (a, b). One therefore can use the coefficients obtained above in the coordinates 
of (e 1 , e 2 ) to obtain the contributions of this ray to the Stokes parameters in the coordinates 
of (a, b). 



2.3. General Relativistic Radiative Transfer 



In magnetized plasmas around black holes, photons experience both the general relativis- 
tic light bending an d plasma birefringence effects and change their propagation directions. 
Fanton et al.l (119971) showed h ow photons propagate along the null geodesies near black holes. 
Broderick fe Blandfordl (120031 ) discussed how the plasma effect causes refraction, which makes 
the photon propagation path deviate from the null geodesies. However, such refraction is 
significant only if either the electron cyclotron frequency z/g or the electron plasma frequency 
up is comparable to the observation frequency z/ b s - In the observational band (from mil- 
limeter to near-infrared band) we are interested in, we always have z^s/Vobs < 10~ 3 , and 
^p/^obs < 10~ 4 . Therefore, although we consider the Faraday conversion and rotation effects 
of magnetized plasmas, we still assume photons of both natural modes propagating along 
the null geodesies. 



We use the ray-tracing code discussed in iHuang et al.l (120071 ) to determine the pho- 
ton trajectory from a specific direction as observed at infinity. Along the trajectory which 
crosses the emission region, we record the four-wave-vector accretion flow velocity u 11 , 
electron temperature T e = ^ c m e c 2 /k 1 number density n, and three magnetic vector B* at 
each line element. We con vert the three magnetic vector B J into four- vector following 
e.g., iGammie et al.l (120031 ) . The four- vectors of the elliptical axe s of synchrotron emis 



sum [e lfJl ,e 2fl ) are calculated the same way as that discussed by iBroderick fe Blandford 



(120041 ) . The four- vectors of the reference coordinates (a M ,6 M ) are calculated according to 



the parallel transport in the general relativistic theory with (a*, b l ) arbitrarily set as (0, 0) 



(1 Chandr asekharl 1 1 9 831 ) . The Stokes parameters are not conserved along the photon trajecto- 
ries due to gravitational effect. However, the photon occupation numbers Ms = (Mj,Mq,Mu,Mv), 
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defined as As = S/u 3 , are Lorentz invariants. Therefore, the radiative transfer equation (J7|) 
can be rewritten as 

-77 — As = £ O/k, (f = 4> £ = *>K), (23) 

d£ i> J 

with the differential affine parameter d£ = ds/^ b s and the emission frequency v = —u^k^. 



SIMULATION RESULTS OF POLARIZATIONS 



In this section, we present our simulation results in detail. Th e data we adopt ar e 



mainly from (mult i- epoch) linear polarizat io n observations r eport ed by 



Bower et al 



Aitken et al. 



Eckart et al. 



( 2000 ): 
(2006). 



20051 ): iMacquart et al.l (120061 ) : iMarrone et al.l (120061 ) . and 
marked with crosses in Panels (a) , (b), a nd ( c) of Figure El Two c ircular polarization data 
in Panel (d) are from iBower et al.l (1200 ll ) and IMarrone et al.l (120061 ) . The ties in X-ray band 



are from iBaganoff et al. 



adopted in iHuang et al. 



1 2001). Other data of luminosity shown in Panel (a) are the same 
fl2008h . 



3.1. Configuration of the Magnetic Field 



We adopt the two-temperature magneto-rotational instability (MRI) driven Keplerian 
accretion flow mo del of magnetized plasmas around the central black hole in Sgr A* [see in 



Liu et al. for details]. A black hole mass of Mbh = 4.1 x 1O 6 M , where M is the 

solar mass , is recently deriv ed from observations of the orbital motion of the short-period 
star SO-2 ( Ghez et al. 20081 ). The corresponding model parameters include the ratio of the 
viscous stress to the magnetic field energy density (3 V [fixed at 0.7, see Pessah et al. (2006)], 
the ratio of the magnetic field energy density to the gas pressure f3 p , the electron heating rate 
indicated by a dimensionless parameter C\, the mass accretion rate M, the inclination angle 
i and the position angle of the axis perpendicular to the equatorial plane of the accretion 
flow. In this model, thermal electrons are energised by sound waves with the heating rate 
given by r^ 1 = Cg/3CiH(v e ), where eg, H, and (v e ) are the sound speed, the scale height of 
the accretion flow, and the mean speed of the electrons, respectively, and we have assumed 
that the scattering mean free path of the electrons by the sound wave s is scaled with H . A 
comparison of this heating mo del with other two-temperature models ( lYuan et al.ll2003l ) can 
be found in iLiu et al.l (j2007bl ) . 



MHD simulations of the MRI show that the magnetic field is dominated by its azimuthal 
component due to the shearing motion of the Keplerian accretion flow. Previous pseudo- 
Newtonian modelling of the sub-millimeter to NIR polarization of Sgr A* has assumed 
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that the magnetic field does not ha ve vertical and radial components (IMelia et al.l 12001 



Bromley et al.ll200ll ; iLiu et al.ll2007al ). Here we consider more realistic configurations. Most 



MRI simul ations adopt poloidal or toroidal magnetic fie ld loops confined to an initial torus 
of matter (IHawley &: Balbud 120021 ; iGammie et al.ll2003l ). Simulations with initial magnetic 
field lines perpendicular to the disk plane tend to produce stronger magnetic fields in the 



final quasi-steady state (IPessah et al. 



20061 ). The field lines are dragged and twisted due to 



the Keplerian motion of the accretion flow. If there is a large scale poloidal magnetic field, 
the field lines will be nearly parallel (anti-parallel) to the velocity field in the upper half to 
the equatorial plane and anti-parallel (parallel) in the lower half. We will adopt such a mean 
magnetic field configuration in the following. 



3.2. The Fiducial Model 



We first obtain a fiducial model to the spectrum, linear polarization (LP) fraction, 
electric vector position angle (EVPA), and circular polarization (CP) fraction with f3 p = 0.4, 



d = 0.47, M = 6xl0 17 g-s 



40°, and = 115°. No external Faraday rotation measure 



is introduced. The corresponding results are shown by the solid lines in Panels (a),(b),(c), 
and (d) of Figure [2], respectively. It is consistent with the fiducial model in Huang et al. 
(2008). This Keplerian accretion flow reproduces the emission in the millimeter to sub- 
millimeter bump of Sgr A*, but underestimates the emission in the NIR and X-ray bands. 
The thick solid lines in Figure [2] (a) represent the synchrotron (to the left) and synchrotron 
self-Comptonization (SSC, in the middle) radiation components, and the thin solid line to the 
right represents the bremsstrahlung radiation. This result is different from that of Liu et al. 
(2007), where they are able to reproduce the observed NIR flux level. In their model, an un- 
polarized j et /outflow comp onent were introduced, which effectively suppresses the LP below 
100 GHz (ILiu et al.ll2007al ). We don't have such a component here. The onset frequency 
of the LP therefore constrains the model parameter space so that the model predicted NIR 
and X-ray fluxes are more than one order of magnitude lower than the observed values. The 
quiescent-state X-ray emission is extended and has been attributed to thermal emission from 
plasmas near the capture radius (IBaganoff et al.ll2003l ; lQuataertll2004l ; IXu et al.l 120061 ) . The 
X-ray emission from our accretion flow at small radii should be lower than the observed 
value. However, there is evidence for a quiescent-state NIR emission from the direction of 
Sgr A* (IDo et al.ll2009l ). This emission can be caused by smaller flares or a quasi- stationary 
emission component. If it is indeed powered by the accretion flow at small radii, it must be 
dominated by emission regions within 6GMbh/c 2 , where G is the gravitational constant and 
c is the speed of light, or associated with outflows, which may also be responsible for the 
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long wavelength radio emissions. @ 

As shown in Figure [2] (b) and (c), the LP fraction is a few percent in the centimeter 
bands. It can be reduced to the observed level either by depolarization of an external medium 
or through the dominance of an un-polarized emission component from a jet /outflow com- 
ponent. In this band, only the red-shifted side of the accretion flow is optically thin. Hence, 
the EVPA is perpendicular to the accretion flow axis projection, since the magnetic field is 
toroidal and the emission is dominated by the extraordinary mode. When the observational 
frequency increases to ~100GHz, the LP degree increases because the Faraday rotation and 
Faraday conversion coefficients decrease with the observation frequency v as v~ 2 and 
respectively. At the same frequency, the near and far sides of the accretion flow become 
optically thin and dominate the polarized emission so that EVPA preforms a ~ 90° flip to 
be parallel to the disk axis projection. In the sub-millimeter band, the LP degree increases 
to ~ 10%, explaining observations well. Furthermore, the EVPA has a small increase of 
~ 20° in this band. With further increasing of the observational frequency, the blue-shifted 
side finally become optically thin at ~ 2 — 3 THz, making the whole disk optically thin 
with its blue-shifted side being the dominant emission region. Thus, the EVPA flips back 
to be perpendicular to the axis projection again. In the quiescent state discussed in this 
paper, there is almost zero polarization in the NIR band because the non-polarized SSC 
emission component dominates. However, the EVPA is still meaningful for flaring activities 
dominated by the synchrotron emission. 

The predicted CP fraction is shown in Figure [2] (d). The CP degree is nearly zero in the 
centimeter bands, so that the observed CP fraction of several percents (Bower et al. 2002) 
might be attributed to a jet/outflow component (Beckert & Falcke 2002). At sub- millimeter 
and shorter wavelengths, however, the CP amplitude increases gradually, then significantly 
to greater than 10% at ~1 THz. We show the CP fraction in different scales for clarity 
of presentation here. Notice that the LP fraction also has some oscillations in the same 
band. These high LP and CP may be due to the combination of the general relativity and 
birefringence effects, and can be tested with future THz polarization observations. 

Interestingly in the millimeter to sub-millimeter band, the CP degree is mostly negative, 
i.e., the polarization is left-handed . Here, we adopt tw o observational data for CPs, a limit 



of -1.8% at 11 2 GHz reported bvlBower et all (120011 ) and a detection of 0.5% ± 0.3% 

at 340GHz by iMarrone et aD toO& \ (both marked with red crosses). Although they men- 



tioned that the detections were quite uncertain due to weather conditions and instrumental 



7 The black hole does not have a spin in our pseudo-Newtonian model. The accretion disk has an inner 
boundary radius of 6GA/bh/c 2 . 
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effects, the left-handed properties seemed to be real. Our model therefore reproduces the 
CP observations. However, we notice that the CP degree is very sensitive to the magnetic 
field structure and model parameters. Thus, our predictions here should be interpreted with 
precautions. 



3.3. The Inclination Angle Dependence 

The dependence of these results on the inclination angle i is also shown in Figure 
[2] with the corresponding models indicated in the legend of Panel (b). The other model 
parameters have been fixed here. Results for i = 20°, 30°, 50°, and 60° are shown with the 
3-dots-dashed, long dashed, dot-dashed, and dotted lines, respectively. Different from the 



pseudo- Newtonian cases [see iLiu et al.l (j2007al )]. the changes in the inclination angle slightly 



affect the centimeter polarization and high frequency flux density. This is due to two effects: 
absorption of emission from the edge of accretion flow and amplification of emission by the 
GR light bending. 

However, the inclination angle affects the LP degree significantly. Generally speaking, 
the onset frequency of high LP decreases and the LP degree in sub-millimeter band increases 
with the increase of the inclination angle. This is because the projected magnetic field lines 
are close to parallel lines in the nearly edge-on case while becoming circular helices in the 
nearly face-on case. Parallel magnetic field lines result in the highest LP degree (up to 30% 
in sub- millimeter band), while circular helices result in a zero LP degree due to the rotation 
symmetry of the flow. The observed LP degree s in the sub-millimeter band remain stable at 



10% since its discovery (lAitken et al.l 120001 ) . which provides a constraint of ~ 40° ± 10 c 



on the inclination angle of the Keplerian accretion flow. 

For the EVPA, it flips at a lower frequency for the first time and a higher frequency 
again. This is because the near side of the inclined disk dominates the emission in a wider 
band. For the i = 60° case, the near side dominates after the first flip, and the second flip 
doesn't appear until the frequency increases to the NIR band. The second flip may disappear 
in the band we are interested in for higher values of the inclination angle. The second flip 
between sub-millimeter band to NIR band is quite important. A stable difference of ~ 80° 
in the EVPAs between 230GHz and A 2.2/im has been observed, which suggests a stable 
geometry for the accretion plasma around Sgr A* if the NIR emission also originates from 
the disk. Comparing the fiducial model and the other four cases shown in Figj2l a mildly- 
inclined i ~ 40° accretion flow appears to better reproduce such a difference i n the EVPA 



caused by changes in optical depth across the emission structure. However, iMeyer et al. 



(|2007l ) derived a highly-inclined i > 70° disk to obtain significantly variable light-curves in 
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the NIR band. The NIR flare emission was associated with a bright spot orbiting the black 
hole at the last stable orbit of a disk with a toroidal magnetic field similar to our magnetic 
field configuration. They also assumed that the sub-millimeter emission comes from a jet to 
explain the EVPA difference, while the jet/outflow component in our model only dominates 
the emission in millimeter and longer wavelengths. Therefore, the difference of inclination 
angle predictions is caused by different assumptions of the geometry for the emitting region, 
which can be checked by future VLBI observations in the sub-millimeter band. 



3.4. The Position Angle Dependence 



The observed EVPA is related to the position angle (G) of the accretion flow axis. 
Observationally, the EVPAs in the sub-millimeter band show frequent variability, presumably 
associated with flares. At 230 GHz, however, it was detected as ~ 1 17° ± 24°, showin g 
relatively stable value s, during several epochs within more than one year (IBower et al.ll2005l ). 
Marrone et al.l (120071 ) also reported comparable values of the EVPA at 230GHz. At A 2.2/xm, 



on th e other hand, the EVPA is reported to be stable with a value of 60° ± 20° (IMeyer et al. 
20071 ). In the fiducial model, we find = 114° so that the predicted EVPA can give a 



best fit to both measurements at 230GHz and A 2.2 /iin. In practice, values in the range 
of ~ 115° ± 20° are all acceptable. Interestingly, we find that the predicted firs t flip in the 
EVPA also fits the m ulti-epoch detections at 86 GHz well (IMacquart et al.ll2006l ). Moreover, 
Marrone et aD ((2003) observed the EVPA at 340GHz simultaneously with that at 230GHz 
and reported an averaged ~ 30° increase from 230 GHz to 340 GHz. Our model also 



reproduces this increase in the sub-millimeter band. IMeyer et al.l (120071 ) derived a O of 
~ 110°, which is consistent with ours, although their emission model is quite different. 

Notice that we have assumed that all the observed EVPAs are intrinsic to the accretion 
flow, i.e., the depolarization and Faraday rotation only occur in the accretion region of Sgr 
A*. It is likely that the emission from the accretion flow may experience depolarization 
and Faraday rotation by an external medium described with a rotation measure (RM) 



4.4 x 10 5 rad • m 2 derived from millimeter to sub-millimeter observations (IMacquart et al. 



20061 ). In Panel (a) of Figure [3], we plot the corresponding EVPAs of three well- fit models 
in Figure [2j with different position angles. Multi-epoch observations of the EVPA from A = 
3 mm to A = 2.2 ^m can be reproduced wi th = 147° for the fiducial model. Fits with 



another RM of 5.6 x 10 5 rad • m 2 derived by iMarrone et al.l (120071 ) from their simultaneous 
observations at 230 GHz and 350 GHz are plotted in Panel (b). = 159° for the fiducial 
model. Obviously, the position angle depends on the external rotation measure. More 
observations in the sub-millimeter and NIR bands will show whether an external rotation 
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measure is necessary. 

3.5. Dependence on the Mass Accretion Rate 

For these radiatively inefficient accretion flows, the gas pressure is dominated by proton 
and ions, which have a relatively high temperature Liu et al. (2007). Pessah et al. (2006) 
showed /3 U ~ 0.7. Then both the viscosity (and therefore the mass accretion rate) and the 
magnetic field energy density are proportional to (3 P . The mass accretion rate therefore 
determines the amplitude of the magnetic field. For a given mass accretion rate, the density 
is inversely proportional to f3 p . Figure H] shows the dependence of the results on the mass 
accretion rate, where (3 P and C\ are also adjusted to reproduce the millimeter and sub- 
millimeter spectrum. For M = 1.5 x 10 18 g • s _1 , (3 P = 0.2, C\ = 0.94. The density profile 
increases by a factor of 5 with respect to the fiducial model. And for a lower accretion rate 
of M = 2 x 10 17 g • s _1 , (3 P = 0.7, C\ = 0.269. The density profile decreases by a factor of 
5.3 with respect to the fiducial model. The inclination angle of the accretion flow is fixed at 
40° and is obtained by fitting observations at 1.3 mm and 2.2/z m wavelengths. It can be 
seen that with the increase of M and therefore the magnetic field B and n, a lower electron 
temperature is needed to reproduce the spectrum. This leads to a lower cutoff frequency 
for both the synchrotron and SSC components. The emission is also strongly self-absorbed 
so that the onset frequency of strong LP occurs at about 150 GHz. Opposite effects can 
be seen with the decrease of M and there is a higher X-ray emission flux due to the SSC. 
However, the X-ray flux is still one order of magnitude lower than the observed quiescent- 
state thermal X-ray emission from a large scale accretion flow (Xu et al. 2005). Even higher 
NIR and X-ray emission may be produced with further decrease of the mass accretion rate 
and increase of the electron temperature. However, radio emission will show strong LP. An 
un-polarized emission component needs to be introduced to suppress the LP to the observed 
level (Liu et al. 2007). 

Moreover, for lower values of M, the high amplitude of the CP degree starts at lower 
frequencies, similar to the high LP degrees, e.g., the dash-dotted line shown in Figure H] (d) 
with (3 P = 0.7. Current observations of the LP and CP therefore constrain the M to be near 
6 x 10 17 g s" 1 with an uncertainty of a factor of ~ 2. 



8 This result is for the pseudo-Newtonian potential without a black hole spin. The range of M will be 
different in a Kerr metric with a significant spin. 
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3.6. Degeneracy Between the Electron Heating and Magnetic Field 

We have shown that the emission spectrum and polarization give good constraints on 
the orientation of the accretion flow and the mass accretion rate. However, the constraint 
on the electron heating rate and the magnetic field parameter j3 p is generally poor even with 
f3 v fixed. Figure shows several reproductions to the observed spectrum and polarization 
with M = 6 x 10 17 g s -1 , p p = 0.2; d = 0.8 (long-dashed lines), and (3 P = 0.7; d = 0.33 
(dashed-dotted lines). The current observations will not be able to distinguish these models. 
However, there are changes in the SSC and bremsstrahlung components. These can be 
understood as the following. 

For a given M, the magnetic field is fixed because both of them are proportional to 
j3 p N Q Tp, where T p is the proton temperature and does not change with (5 P . With the increase 
(3 P , we increase the viscosity and therefore the radial velocity. The density will decrease as 
N oc l//3 p . For a given optically thin synchrotron luminosity required to reproduce the 
millimeter and sub-millimeter spectral hump, N Q B 2 T 2 should not change. Therefore the 
electron temperature T e should be proportional to f3 p . Then with the increase of (3 P , the 
synchrotron spectrum should be slightly harder, i.e., its spectral cutoff oc BT% oc (3 P should 
shift toward higher frequencies. We therefore expect less low frequency SSC flux density due 
to the decrease of the electron density. However the SSC component should cut off at a higher 
frequency, which is proportional to T 2 oc (3 p . The SSC luminosity, which is proportional to 
NqT 2 , however does not change. To make the SSC cut off in the Chandra X-ray band, we 
will need to increase (3 P by a factor of ~ 10. However, this will reduce the flux density 
by the same factor. Therefore the SSC component cannot have significant contribution to 
the observed quiescent state X-ray emission. The bremsstrahlung flux density scales with 
NZ/Tl' 2 oc (3 P 9/4 . A very low value of (3 p is required for this component to have significant 
contribution in the X-ray band. There is therefore a degeneracy between the electron heating 
rate and the magnetic parameter. By introducing a spin parameter to the black hole, the 
emission volume will be reduced and therefore enhancing the SSC X-ray flux. One may be 
able to break this degeneracy with the quiescent state NIR and X-ray flux densities. 



3.7. Polarized Images in the Millimeter and Sub-Millimeter Bands 

Sgr A* is embedded in the inter-stellar medium (ISM) on the Galactic plane. Emission 
from the accretion flow will experience significant scattering, which smoothes the image 
without changing its total flux density and polarization properties. We adopt an elliptical 
Gaussian structure for the scattering screen with a full width of half maximum (FWHM) 
along its major and minor axis being $ ma $ = (1.39 ± 0.02)(A/1 cm) 2 mas and $ m i n = (0.69 ± 
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0.06) (A/1 cm) 2 mas, respectively, and a position angle ~ 80° (IShen et all 120051 ) . The size of 
the scattering screen becomes smaller with the decrease of the observation wavelength. In 
practice, the black hole shadow structure is significantly washed out in the millimeter and 
longer wavelength band by this scattering. However, the scattering becomes less important 
in sub-millimeter band. We plot images of the total emission I, LP :v /Q 2 + [/ 2 , and CP: 
V of the fiducial model at several wavelengths A = 3.5mm, 1.3mm, 0.86mm, and 0.6mm in 
Figure [61 The black dashes in the LP emission represent the EVPAs % The left-handed 
and right-handed regions in the CP emission are respectively represented with grey and red 
color. 

The image of the total emission at A = 3.5 mm is smeared by the screening and has an 
~ 80°-oriented elliptical Gaussian structure with the major FWHM of ~ (0.189±0.003)mas. 
T his size is consisten t with but somewhat smaller than the observation of ~ 0.21 mas reported 



m 



Shen et all (120051 ). The LP and CP emission at A 3.5 mm are also smoothed significantly. 
The spurs mark the EVPAs, mostly in the radial direction as a result of the toroidal magnetic 
field configuration. The CP emission is dominated by a grey region representing a left-handed 
polarization. 

Clear black hole shadow structures can be found in images of the total emission at 1.3 
mm or shorter wavelengths. Interestingly, clumpy patterns are seen in the sub-millimeter 
LP with bright and faint regions, and in the CP emission with their left-handed and right- 
handed regions. Such complex patterns are results of the gravitational light bending and 
the birefringence effects posing new challenges for future polarization-sensitive VLBI obser- 
vations. 



4. SUMMARY AND DISCUSSION 

In this paper we explore the parameter space of the two-temperature MRI driven Kep- 
lerian accretion flow model in a non-spin pseudo-Newtonian potential for the time averaged 
millimeter and shorter wavelength emission from Sgr A* (Liu et al. 2007). The model repro- 
duces the observed emission spectrum and polarization with an inclination angle of 40 ± 10° 
and a mass accretion rate of ~ 6 x 10 17 g s -1 . The former is mostly determined by the 
amplitude of the LP fraction and the latter is well-constrained by the onset frequency of 
prominent LP. The LP is low in the centimeter wavelength due to strong self-absorption in 
the emission region. The orientation of the accretion flow projected on the sky depends on 
the amplitude of an external Faraday rotation measure. And the model predicted NIR and 
X-ray fluxes are more than one order of magnitude lower than the observed low flux levels. 
Although nearly identical millimeter and sub-millimeter spectrum and polarization can be 
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obtained by adjusting the electron heating rate and the magnetic parameter, none of these 
models can produce NIR and X-ray fluxes comparable to the observed values. This strongly 
suggests that the black hole is rotating so that the last stable orbit can be smaller, resulting 
in more SSC and bremsstrahlung emission. We are in the process of this study and the 
results will be reported in a separate paper. 

The time averaged source size measurements provide critical constraints on the model. 
Although we showed that our model are consistent with current observations by simulating 
the images of the accretion flow at different wavelengths. Given the challenges in imaging 
the black hole with the global VL BA, a direct fit to the observed visibility is needed to fully 



utilize the observed information (IDoeleman et al.ll2008l ). 



Very rich information is contained in the flare observations. To fully explore the impli- 
cation of these observations, one needs to carry out GRMHD simulations of the accretion 
flow. Our currently study of the time averaged properties will lead to good constraints on 
the properties of the plasmas near the black hole. These will be helpful for setting up the 
simulations to recover the essential observations. Besides GRMHD simulations, the kinetics 
of electron heating and acceleration must also be addressed (Liu et al. 2004). The GR 
ray-tracing treatment of polarized radiation transfer we have here in this paper, GRMHD 
simulations, and a kinetic model for the electron acceleration should be combined to develop 
a self-consistent model for the emission structure near the black hole. In light of continuous 
high resolution and sensitivity observations, these theoretical developments will be essential 
to uncover the nature of the black hole and its interaction with the plasmas feeding it. 
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Fig. 1. — Sketch map of polarization in synchrotron radiation: vectors k, B, e 1 , and e 2 
represent the wave vector, magnetic field, e-mode, and o-mode, respectively. \ an d d\ 
respectively represent the EVPA defined in arbitrary coordinate system (a, b) and the change 
of EVPA due to radiation transfer effect. See text for detail. 
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Fig. 2.- 
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( 1200ll ); Marrone et al.l (120061 ). iBaganoff et al.l (120011 ) .etc. . (a) spectrum of synchrotron, SSC 
radiation (thick) and bremsstrahlung radiation (thin); (b) linear polarization degrees; (c) 
EVPA; (d) circular polarization degrees. 
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Fig. 3. — Results of EVPA three mode ls in Figj2] with external rot ation measure of (a) 
4.4 x 10 5 rad • m" 2 flMacquart et alJbood ) and (b) 5.6 x 10 5 rad • m" 2 jMarrone et al.l l2007r > 
considered. 
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Fig. 4. — Simulated results for models with (3 P = 0.2, M — 1.5 x 10 18 g s 1 (long-dashed), 
P P = 0.4 (solid), and f3 p = 0.7, M = 2 x 10 17 g s" 1 (dash-dotted), where fi p -C x = 0.4 • 0.47, 
and the luminosity in sub-millimeter band does not change significantly. 
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Fig. 6. — Images of total intensity (J), LP emission (^/O^PU^), and CP emission (V) 
predicted for Sgr A* at observational wavelengths of 3.5 mm, 1.3 mm, 0.86 mm, and 0.6 
mm, from top to bottom, respectively. In images of LP emission, black spurs represent 
the averaged EVPAs in the nearby regions. In images of CP emission, the left-handed and 
right-handed regions are respectively presented in grey and red. 



